Optimal coherent control of CARS: signal enhancement and background elimination 
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The abihty to enhance resonant signals and eliminate the non-resonant background is analyzed for 
Coherent Anti-Stokes Raman Scattering (CARS). The analysis is done at a specific frequency as well 
as for broadband excitation using femtosecond pulse-shaping techniques. An appropriate objective 
functional is employed to balance resonant signal enhancement against non-resonant background 
suppression. Optimal enhancement of the signal and minimization of the background can be achieved 
by shaping the probe pulse alone while keeping the pump and Stokes pulses unshaped. In some 
cases analytical forms for the probe pulse can be found, and numerical simulations are carried out 
for other circumstances. It is found that a good approximate optimal solution for resonant signal 
enhancement in two-pulse CARS is a superposition of linear and arctangent type phases for the 
pump. The well-known probe delay method is shown to be a quasi-optimal scheme for broadband 
background suppression. The results should provide a basis to improve the performance of CARS 
spectroscopy and microscopy. 

PACS numbers: 



I. INTRODUCTION 



Coherent anti-Stokes Raman scattering(CARS), as a 
four- wave nonlinear process^, has been widely used in 
the past few decades to study chemical systems in so- 
lutions, reactions in the gas phase, and vibrational dy- 
namics in gas and condensed phases. CARS microscopy 
is a recently implemented technique for imaging biologi- 
cal species, which was pioneered by Duncan et al. using 
two dye lasers and developed^ by Xie et al. for high- 
sensitivity applications^'^. As a combination of ultra- 
fast nonlinear spectroscopy and microscopy, CARS mi- 
croscopy is a highly chemically selective and sensitive 
technique that employs a CARS signal of an unlabeled 
sample and provides higher spatial resolution than two- 
photon fluorescence microscopy. 

When imaging biological samples, the typical width of 
Raman transitions is a few wave numbers, so picosec- 
ond pulses are widely used in CARS spectroscopy and 
microscopy^' Not only is the narrow bandwidth of 
the picosecond pulse adequate to detect specific Raman 
bands, it also produces low nonresonant background. On 
the other hand, when investigating broadband CARS 
spectra^'^, especially in the Raman fingerprint region 
which spans from 800 to 1800 cm~^, it is necessary to 
bring in broadband femtosecond pulses. Femtosecond 
CARS can be employed not only for direct imaging but 
also as a tool to determine some microscopic and macro- 



scopic parameters, such as molecular anharmonicity^^ 
and temperature^^'^^. However, this situation creats a 
dilemma when the pulse bandwidth is ~ 1000 cm~^, as 
the nonresonant background becomes significant to possi- 
bly submerge the resonant signal and the fine vibrational 
structure of CARS^^'^^. The large nonresonant contribu- 
tion affects the shape of CARS spectra and complicates 
data analysis, which becomes an obstacle for femtosec- 
ond CARS. Hence the study of signal enhancement and 
background suppression is important for effective Ra- 
man selective excitation in CARS, in which a broadband 
femtosecond pulse is employed to excite multiple Raman 
modes^^^^. Resonant signal enhancement and nonreso- 
nant background suppression of CARS have been studied 
for many years. Polarization CARS adjusts the polariza- 
tion of the pump and Stokes pulses to suppress nonreso- 
nant background"*^^^"*^. Time-resolved CARS^^ uses tem- 
porally overlapped pump and Stokes pulses along with 
a delayed probe pulse to generate a signal. This pro- 
cedure can also eliminate the nonresonant background, 
which is essentially instantaneous while the resonant sig- 
nal has a much longer decay time. Scully et al. pro- 
posed hybrid CARS^^, in which the broadband pump and 
Stokes pulses produce maximal Raman coherence and 
the narrow-band time-delayed probe pulse suppresses 
the nonresonant background. Cheng et al. reported 
that epidetected CARS microscopy can significantly re- 
duce the solvent background^. The advent of spatial 
light modulators (SLMs) has enabled coherent control 
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and mode-selective excitation of CARS with femtosecond 
pulses^^^''. Background suppression without loss of the 
resonant signal has also been explored widely with the 
development of femtosecond phase shaping techniques. 
By using SLM phase modulation and a variable wave 
plate, Silberberg et al. combined phase and polarization 
control to yield background-free single-pulse multiplex 
CARS spectra with high spectral resolution^^. Leone 
et al. combined interferometric and polarization/phase 
control to demonstrate a method of single pulse inter- 
ferometric CARS spectroscopy^^'^^, which could extract 
the imaginary and real parts of the background-free res- 
onant CARS spectrum in a single experimental measure- 
ment. Offerhaus et al. employed vibrational phase con- 
trast CARS, in which the measured phase components 
in the focal volume allows enhanced sensitivity and in- 
creased selectivity^^'^^. With the optical fields driving 
the CARS process and the local oscillator used for het- 
erodyning both derived from a single beam by pulse shap- 
ing, Motzkus et al. proposed highly sensitive single-beam 
heterodyne CARS microspectroscopy^^, in which the sen- 
sitivity of chemically selective detection at microscopic 
resolution is dramatically increased. It has also to be 
mentioned here that while phase shaping is effective to 
suppress nonresonant background in broadband CARS, 
there are also post-processing approaches proposed by 
Bonn^^^^ and Cicerone^^'^^ to retrieve background-free 
and noise-reduced CARS spectra. 

In this work, we explore the optimal control of the sig- 
nal and background of CARS with various phase shaping 
schemes. The paper is organized as follows: Section II 
shows the optimal phase shaping schemes at a specific 
frequency. Section III investigates the control strategies 
for broadband CARS and the conclusions are given in 
Section IV. 



II. LOCAL OPTIMAL CONTROL 



CARS is a four-wave nonlinear process as shown in 
Fig.(l). Three laser pulses are used to produce the 
CARS signal: the pump pulse Ep{ujp)^ the Stokes pulse 
Es{ujs)^ and the probe pulse Epr{ujpr). The CARS signal 
IcARS as) is a coherent superposition of resonant third 

) and nonresonant 



(3) 

order nonlinear polarization Pr (cj, 



third order nonlinear polarization Pnr (<^, 
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and 



X Ep{ujp)E*{u;s)Epr{Cjpr)S{Cjas—ojp-\-u;s—oJpr) 



XEp{0Jp)E* {Cjs)Epr{0Jpr)5{0Jas—0Jp+0Js—UJpr) 



(1) 



(2) 



(3) 



here is the Raman frequency between energy level 
1 and 2, 2r is the level width, C is a constant which 
depends on the material property, and Xnr is the non- 
resonant third-order susceptibility. The resonant sig- 
nal intensity at the frequency uJas is \Pr^\oJas)\'^ ^ the 
nonresonant background intensity is |Pn?(^as)P, and 
the integrated CARS intensity is / = /q^^ \Pr^\oJas) + 

Pr^^\uJas)\^dCUas^ 

In our work, the carrier frequencies of the pump. Stokes 
and probe pulses are denoted as Qp, fts and flpr, re- 
spectively. Then if all the pulses are transform limited 
pulses(TLPs), the peak of the resonant and nonresonant 
signals in the spectrum are both located at the frequency 

^as — ^p 

A Gaussian amplitude profile in the frequency domain 
is used for the pump. Stokes and probe pulses in our 
theoretical treatment and simulations. 



Ek (a)fc)= i 



fc-f^fc), k={P, S, Pr}, (4) 
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Figure 1: Energy level diagram of the CARS process. The left 
panel corresponds to the resonant signal generation: pump 
and Stokes pulses generate coherence between two vibrational 
levels, when they have a frequency difference which coincides 
with the Raman resonance Qr. The probe pulse then induces 
the anti-Stokes signal. The right panel corresponds to the 
nonresonant background contribution: the nonresonant back- 
ground is produced via an intermediate virtual state that does 
not reflect the resonant molecular energy level. 



where 2vlri2A/e and 2vln2A/e are the corresponding 

spectral full widths at half-maximum (FWHM), and 

^k{^k — ^k) is the frequency-domain phase profile. For 

simplicity, a frequency variable translation, uj^ = ^k—^k^ 

is made {i.e. ^k{^k) = ^k{(^k - ^k) and Ek{oJk) = 

Ek{CJk) = ^e-^ll^le'^^^^^^), and P^^\^as) is writ- 
^k 

ten as P^. 

There are several experimental configurations for 
CARS: single-pulse CARS^ (all three photons required 
are supplied by the same short optical pulse), two-pulse 
CARS (the pump pulse also acts as probe pulse in the 
experiment) and three-pulse CARS. In this section, opti- 
mal control at a single signal frequency is investigated in 
two-pulse and three-pulse CARS. Optimal control strate- 
gies with the probe pulse are also explored. Only the in- 
tensity at a specific frequency (Co as = ^as) is considered 
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throughout this section, and broadband optimal control 
will be discussed in Section III. 

For simplicity, unless otherwise stated, Ap = A^^ = 
A5 = A is assumed in the following without loss of gen- 
erality, as the conclusions still hold when the bandwidths 
are different. 



Control of the signal or background by shaping 
the probe pulse in three-pulse CARS 
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Figure 2: The resonant signal and nonresonant background 
with different phase shaping schemes for the probe pulse while 
keeping the pump and Stokes pulses unshaped in three-pulse 
CARS. The arctan(c(;pr/r) phase (red solid lines) , vr-step 
phase (blue dash dotted lines) and TLP (black dashed lines) 
shaping schemes are shown together for comparison. The 
bandwidths of the pump, Stokes and probe pulses are the 
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: 50cm 



Optimal control for signal enhancement or background 
suppression with the probe pulse has been discussed 
previously^^^^, so only a brief description is given here. 
Fig. (2) shows the optimal control strategies: when the 
phase of the probe pulse is either the arctan(c(jp^/r) or 
a TT step phase function, the resonant and non-resonant 
signals achieve their maximal and minimal values, respec- 
tively. The optimality condition can be gained analyti- 
cally as following: 

Eq.(2) leads to 



(5) 



: exp 



where a{ujpr) = — arctan(cjpr/r) + 7r/2 is confined in 
the domain [0,7r]. It is easy to see that |P^p reaches 
its maximal value when the phase contribution to the 
integrand exp [i^pr{oJpr) + ia{ujpr)] becomes a constant, 
i.e. 



^pr{(jJpr) = arctan(a;pr/r) + constant, 
with the resultant maximal peak intensity 



(6) 



(7) 



Here K is a modified Bessel function. The condition 
for maximal or minimal resonant signal intensity could 
also be established by using the variational method (see 
the appendix A), which may aid in exploring the control 
landscape^^ for CARS. 

The non-resonant background can be derived from Eq. 
(3) 



exp[i^priu;pr)] (8) 



Thus |Pnr(^as)P rcachcs its minimal value of zero 
when ex.-p[i^pr{oJpr)] is an anti-symmetric function, e.g. 
^pr{oJpr) is the TT stcp phasc function about ujpr = 
r 3uj^ 1 

(note that exp — ^ symmetric positive-definite 

function). 

As seen in the bottom panel of Fig. (2), the tt step 
phase function can only eliminate the local component of 
background around Uas = ^as to form a dip in the spec- 
trum. Thus the phase function obtained by minimizing 
\Pnr{^as)\'^ IS locally Optimal, but the background may 
still affect the resonant signal away from the position of 
Coas = ^as in the spectrum. Hence, a broadband back- 
ground suppression method is necessary for CARS, which 
will be discussed in Section HL 



B. Control in two-pulse CARS 

For practical considerations, in many CARS experi- 
ments, the pump pulse also operates as the probe pulse. 
In this subsection, we will concentrate on two-pulse 
CARS, in which only the pump pulse is phase shaped. 
From Eq. (2), it follows that 

P^^\nas) = S-aa cop^-iT ^pi^pr-) [j-oa Ep(uJp) (uJp-\-UJpr)duJp]duJpr- 



2(u;p+a;p^/2)^ 

a"2 



doup 



dujpr- 



(9) 



It is difficult to obtain an analytic optimal phase 
function for \P^{Qas)\' Instead, a numerical solution 
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is presented in Fig. (3), and it shows that the opti- 
mal phase function is approximately a superposition 
of the linear and arctan(cc;p/r)/2 phases: the optimal 
phase is quasi-linear away from ujp = and simi- 
lar to arctan(a;p/r)/2 around ujp = 0. This behav- 
ior may be understood as follows: since the zero phase 
profile (a special case of linear phase) is optimal for 
I ^_^(-p+-pr-/2) +i<i>^(u;^)^^^| ^ arctan(cjp/r) 

phase is optimal for | ^ '^^^^''^^'^^ dujprl , 

then an approximate superposition of the linear and 
arctan(a;p/r) phase functions forms the optimal solution 

for \Pr {ftas)\j which is the combination of these two 
integrals. With no explicit analytical solutions, optimal 
control in two-pulse CARS has to be studied numerically 
case by case. In an experiment, an approximate super- 
position of a step and a proper time-delayed phase could 
effectively produce a large resonant signal^^ in the case 
of small r. 
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Figure 3: Optimal phase function(via the CMA-ES optimiza- 
tion method^^'^"^) for the pump pulse in two-pulse CARS. The 
parameters: Ap = = 50cm~^, F = 4.8cm~^. In the top 
panel, the red solid line is the optimal phase function of the 
pump pulse for maximal resonant signal intensity, the ma- 
genta dotted line corresponds to the arctan(c(;p/F)/2 phase, 
and the blue dashed line is a linear phase profile. The bottom 
panel shows the outcome of optimal pulse (red solid line) and 
TLP (black dashed line) in the time domain. 



To check the effect of this optimal phase scheme in two- 
pulse CARS, the resonant signal and nonresonant back- 
ground spectra with other shaping schemes are shown 
together in Fig. 4 for comparison. As can be seen, \Pr\ 
for the optimal scheme looks similar with that for the 
arctan(cjp/F) scheme. The resonant signal decreases a lot 
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Figure 4: The resonant signal and nonresonant background 
spectra with different phase shaping schemes for the pump 
pulse in two-pulse CARS. The optimal phase (red solid lines), 
arctan(cL;p/r) phase(blue dotted lines), tt step phase (magenta 
dash dotted lines) and TLP (black dashed lines) schemes are 
shown together for comparison. The parameters are the same 
as in Fig. 3. 



with the TT step phase. Compared with the TLP scheme, 
all the other three schemes can suppress the background, 
but only the optimal phase scheme enhances the resonant 
signal. This is easy to understand since our optimized 
objective is the resonant peak signal at uoas = 0. 



C. Control in three-pulse CARS 

In Sec. 11(A), it was shown that the phase function 
arctan(a;p^/r) generates a maximal resonant signal in- 
tensity when only the probe pulse is shaped. In three- 
pulse CARS, all the three pulses can be shaped. Thus 
in this subsection, we will demonstrate if the three-pulse 
shaping scheme can achieve better performance than the 
probe-only shaping scheme in three-pulse CARS. The an- 
alytical and numerical results will show that the config- 
uration of TLPs for the pump and Stoke pulses and an 
arctan(a;pr/r) phase profile for the probe pulse is optimal 
to maximize the resonant signal. The CARS spectra of 
\Pr\ and \Pnr\ with this optimal phase scheme and other 
schemes can be found in Fig. 2. 

From Eq. (2), we have 



(10) 
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where the Raman excitation term is 



A{uJp^UJpr) - 

It is easy to see that 

\A{uJp,UJpr)\ < 



-^2 ^'i^pi^p)-'i^siu;p+ujpr^) 

(11) 



p+u;p^/2)^ 





duJn 



A (12) 



The equahty holds only when ^p{ujp) — ^s{^p-\-^pr) does 
not depend on variable ujp for arbitrary ujpr. There are 
only two cases satisfying this condition: 1) ^p{ujp) and 
^s{^p -\-^pr) are both constant functions^ i.e. the pump 
and Stokes pulse are TLPs; 2) ^p and are both linear 
functions with the same slope: ^p{ujp) — ^s{^p + ^pr) = 
kujpr-\- constant, which is equivalent to linear phase shap- 
ing (or a time delay) scheme for the probe pulse. The first 
case is just the second one with zero slope (no delay). 
Thus, the following equation can be derived 



^3/2 J-oo ujpr-ir 



3u. 



A[UJp,UJpr)aU!pj 



\\A{uJp,UJpr^)\duJpr 



^3/2 J-oo ^ oJp^ — zT I V 2 P 

(when ^p and are constant) 



<c 



e%"BesselK(^0, 



3r^ 

4A2 



(when ^pr{ujpr) = arctan(a;pr/r)) 



(13) 



and the equality only holds when ^p and are constant 
and ^pr{uJpr) = arctan(cjp^/r). When Ap, Apr and A5 
are different, this conclusion still holds. Hence the maxi- 
mal resonant signal is only achieved when the pump and 
Stoke pulse are unshaped TLPs, and the phase of the 
probe pulse takes the form ^pr{oJpr) = arctan(cc;p^/r). 
This fact means that a two-pulse or three-pulse shaping 
scheme is not necessary to achieve an optimal resonant 
signal, which was also verified by numerical simulations 
as shown in Fig. 5. The CMA-ES algorithm, which is 
reliable for global optimization, is employed in the nu- 
merical optimization. In agreement with the above ana- 
lytical result, it is found that the configuration of TLPs 
for the pump and Stoke pulses and the arctan(cjp^/r) 
phase profile for the probe pulse is optimal to maximize 
the resonant signal. 



D. Control of the signal-background difference by 
shaping the probe pulse 

In the subsections above, the optimization of the reso- 
nant signal and non-resonant background is treated sep- 
arately. In the laboratory, however, the signal and back- 
ground are always detected together. As they can not 
simultaneously reach extreme values, it is necessary to 
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Figure 5: Numerical optimal phase functions for achieving 
a maximal resonant signal with three-pulse CARS using dif- 



ferent pulse bandwidths: Ap — 125cm 
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The top and bottom panels correspond to the 



frequency and time domains, respectively. The arctan(c(;pr/r) 
phase (black solid lines) is also shown for comparison. 



study the balance between resonant signal enhancement 
and non-resonant background suppression, which is a 
multi-objective optimization problem. In this subsec- 
tion, we will show how to achieve this goal by shaping 
the probe pulse. 

For this multi-objective problem, it is natural to con- 
sider the optimization of the signal-to-background ratio. 
However, the ratio {\Pr\/\Pnr\) is not a good choice: it 
could become infinite when \Pnr\ = 0, no matter how 
small \Pr \ is. In this work, the difference of the resonant 
signal and non-resonant background intensities is chosen 
as the objective functional. 



J 



k\Pn 



(14) 



where k is the weight factor. This objective functional 
J is the balance of the minimization of |PnrP and maxi- 
mization of |P^p. By maximizing the intensity difference, 
a large signal-to-background ratio can be achieved with 
significant resonant signal intensity. 

According to the variational method, the necessary 
condition for a stationary point of J is 



SJ 



6^pr{uJpr) 



Numerical 



0. (15) 

results indicate that the optimal phase 
^pr{ujpr) is an odd function of uOpr- After a detailed 
analysis shown in appendix B, it is found that the op- 
timal phase for maximizing J is a modified arctan-type 
function. 
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Figure 6: The top panel: Numerical optimal phase function 
of the probe pulse for \Pr\'^ — k |Pnr|^ with different weights 
k. The color of the lines indicates the value of k, which is 
represented in the color bar on the right corresponding to 
\ogiQ{k + 0.1). All the phase functions in this figure could 
significantly suppress the background. The bottom panel: 
The Pareto surface for the optimization of signal enhancement 
and background suppression. With different weights k, the 
value of \Pr\^ is bounded in [0.765, 0.828], while \Pnr\^ is 
always much smaller than |Pr|^. 



where A = k{xnr/Cy — ftpr^ 7 (dependent 

on A) is a parameter determined in Eq. (B9) of appendix 
B, and 6> is a trivial phase angle. When the weight factor 
/c = 0, the optimization of J reduces to the maximization 
of resonant signal, and its solution arctan(cjpr/r) is also 
consistent with the result in Eq. (6). With parameters 
C = 1 and Xnr = 0.1, the numerical optimal phases for 
maximizing |P^|^ — /c |Pnr|^ with different k are shown in 
the top panel of Fig. (6). It is easy to see that the Pareto 
surface of J in the bottom panel of Fig. (6), which is de- 
fined as the set of optimal points (|P^(/c)|, |Pnr(^)P) for 
maximal J with different /c, has two limit points. In the 
first limit, the pure maximization of |PrP without con- 
sidering |PnrP leads to |P^p reaching its maximal value 
of 0.828 while |PnrP has a considerable value of 0.19, 
and the phase function is arctan(cjpr/r)- If Xnr is large 
in this case, |PnrP is large as well, e.g. when Xnr = 1-0, 
then |PnrP = 1-9. Hence the choice of k = will not 
achieve an optimal balance of resonant signal enhance- 
ment and non-resonant background suppression in gen- 
eral cases with a large background. In the second case 



when k becomes large, then |P^p approaches its lower 
limiting value of 0.765 while |PnrP is very small, and 
the phase function also approaches a limiting curve as 
indicated with the red color in the top panel of Fig. (6). 
In this case, the background could almost be eliminated 
while keeping a considerable resonant signal intensity. 
Sec. 11(A) showed the tt step phase could also elimi- 
nate the background, but it produced a relatively small 
resonant signal intensity |PrP = 0.45(with the same pa- 
rameters as here). Hence, the tt step phase is not a good 
Pareto optimal solution for the optimization of J, al- 
though it is widely used in experiments. 

From the analysis above, a practical choice for k can be 
made: 1) when Xnr/C is large {i.e. a large background), 
it is better to choose a large value of /c, which will give a 
solution close to the limit point on the right side of the 
Pareto surface; 2) when Xnr/C is small, a small value of 
k is appropriate. For example, when Xnr/C approaches 
zero, the k=0 is the best choice, in which case there is no 
need to consider further reducing background. 



III. BROADBAND OPTIMAL CONTROL 

In Section II, resonant signal enhancement and non- 
resonant background suppression at a specific frequency 
was studied. However, in many experiments, a large res- 
onant signal and low non-resonant background over the 
entire spectrum is designed, so the local optimal control is 
not adequate. In this section, we will show how to achieve 
broadband optimal control. The tt step phase, which 
achieves perfect local elimination of the background, will 
be investigated again, and a multi-7r step phase scheme is 
further proposed to yield better broadband background 
suppression. Finally, a global numerical optimized phase 
profile for broadband background suppression will be 
shown by employing a specific objective functional. 



A. Multi-TT step phase scheme 

It is shown above that the tt step phase (Heaviside(co'pr)- 
tt) could completely eliminate the background at 
l^as^which can also be achieved at another specific fre- 
quency by changing the tt phase step position. The ef- 
fects of the TT step phase function are shown in Fig. (7). 
As can be seen, this phase scheme could not effectively 
eliminate the background at other frequency components 
at the same time, although it produces a good signal-to- 
background ratio at ^as^ which means that the tt step 
phase is not a good choice for broadband background 
elimination. 

Since a single tt step phase could eliminate the back- 
ground locally and its effective region for background 
elimination is only about A/2 as shown Fig. (7), it is 
natural to construct a multi-7r step phase with jumps 
at several positions for broadband background suppres- 
sion. Every tt phase step corresponds to a local region 
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Figure 7: The top panel: The resonant signal(red dashed 
Une), background (blue dotted line), and the whole CARS 
signal (black solid line) with the tt step phase scheme; The 
bottom panel: The tt step phase profile of the probe pulse 
which steps about uopr — 0. Parameters: A = 50cm ~^ , 
r = 4.8cm-\ 



Figure 8: The top panel: The resonant signal intensity \Pr\ 
and non-resonant background \Pnr \ intensity with the multi-7r 
step phase scheme (red lines) and the time delay scheme (blue 
lines). The bottom panel: The multi-vr step phase profile (red 
solid line) and the time delay phase profile (blue dashed line) . 



for background suppression. Fig. (8) illustrates a multi-7r 
step phase scheme which consists of eight tt steps with an 
equivalent spacing of A/2. With this ladder-like multi-7r 
step phase, the background is perfectly suppressed and 
the resonant signal remains large. The ladder shape of 
the multi-TT step phase is close to a time delay scheme as 
shown by the two curves in the bottom panel of Fig. (8). 
In fact, similar CARS spectra are yielded by these two 
schemes as shown in the top panel of Fig. (8), which 
indicates that the time delay scheme may be good for 
broadband background elimination. This will be further 
demonstrated in the following subsection by numerical 
global optimization. 



B. Numerical global optimization 



Since it is difficult to obtain an analytical optimal 
phase function for broadband background suppression, 
a global numerical optimization of the phase profile with 
a specific objective functional is performed, which is de- 
fined as the difference of the integrated resonant and non- 
resonant signal intensities. 



r-\-oo 

Ir = 



(17) 



/ \Pr{Cjas)?dLOas, (18) 

— oo 

/ + 00 
\Pnr{^as)?d(jJas, (19) 
-oo 



The optimization of this objective functional aims for a 
large resonant signal while suppressing the background. 
Numerical optimization with respect to all three pulses 
is carried out with A = 50cm~^ and F = 4.8cm~^. 
The optimal phase shaping configuration consists of un- 
shaped pump and Stokes pulses and a shaped probe pulse 
as shown in Fig. (9). It can be seen that the optimal 
phase profile is quasi-linear. In other words, the time 
delay scheme, a well-known method for the background 
suppression, is quasi-optimal for broadband background 
elimination with the objective functional defined in Eq. 
(17). 

To check exactly how much improvement can be 
made using the optimal phase scheme, the corresponding 
CARS spectra of resonant signal |P^| and nonresonant 
background \Pnr\ with different phase shaping schemes 
are shown together in Fig. (10). Compared with the TLP 
phase scheme, both the optimal and time delay schemes 
can suppress the background to a low level across a broad 
frequency band and enhance the resonant signal around 
uOas = 0. The spectra of |P^| with the optimal and time 
delay phase schemes are almost the same except the tails 
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away from uOas — 0, while the spectra of \Pnr\ show 
that the background is suppressed more with the opti- 
mal phase scheme especially around ujas = 0- However, 
this difference in the suppression of \Pnr \ is small, which 
leads again to the conclusion that the time delay scheme 
is quasi-optimal. For other types of objective functionals 
which place more importance on broadband background 
elimination, this conclusion still holds. Although the 
time delay scheme could not yield a narrow-band spec- 
trum, it has almost the best performance in suppressing 
the broadband non-resonant background, which is espe- 
cially advantageous when the background is very large. 
In practice, appropriate combinations of the time delay 
and other phase shaping schemes could achieve enhanced 
performance for signal-background control in CARS. 




Figure 9: Optimal phase profile (with the CMA-ES algo- 
rithm) for the maximization of Ir — Inr • Because the numer- 
ical optimal phase function for the pump and Stokes pulses 
are zero phase, i.e. TLP, they are not shown in this fig- 
ure. The top and bottom panels show the phase function 
and amplitude for the probe pulse (red solid lines) in the fre- 
quency and time domain, respectively. The unshaped TLP 
(blue dashed line) and quasi-optimal time delay (black dash 
dotted lines) schemes are also shown for comparison. Param- 
eters: A — 50cm~^ and F — 4.8cm~^ . 



IV. CONCLUSIONS 

Detailed investigations of coherent control for resonant 
signal enhancement and non-resonant background elim- 
ination of CARS via phase shaping schemes lead to the 
conclusion that the maximal resonant signal and minimal 
background at a specific frequency may be achieved by 
shaping the probe pulse only, while keeping pump and 
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Figure 10: The resonant signal and nonresonant background 
with different phase shaping schemes for the probe pulse, 
while keeping pump and Stokes pulses unshaped. The op- 
timal shaping scheme for maximal Ir — Inr (red solid lines), 
quasi-optimal time delay scheme (black dash dotted lines) and 
TLP scheme (blue dashed lines) are shown together for com- 
parison. Parameters are the same as in Fig. 9. 



Stokes pulses unshaped. The optimal probe phase func- 
tion in two-pulse CARS is approximately a superposition 
of linear and arctangent type phases for the pump, which 
enhances the resonant signal more than other schemes. 
As a balance of resonant signal enhancement and non- 
resonant background suppression, the optimization of the 
objective functional |Pr(^as)P — ^\Pnr(S^as)^ could si- 
multaneously generate a CARS signal with large resonant 
component and small background. To achieve broad- 
band non-resonant background suppression, the differ- 
ence of the integrated signal-background intensity over 
the entire spectrum is taken as the objective functional. 
It is found that the optimal phase shaping configura- 
tion consists of unshaped pump and Stokes pulses and 
a quasi-time-delay probe pulse. Numerical simulations 
show that the background is suppressed more especially 
around uJas = with the optimal phase scheme than with 
the time-delay scheme. But the difference in the sup- 
pression of \Pnr\ is small, which leads again to the con- 
clusion that the well-known time delay scheme is a good 
approximation for optimal resonant signal enhancement 
and broadband background suppression. It is expected 
that performing coherent control of the resonant signal 
and background could help improve the performance of 
CARS spectroscopy and microscopy, especially when us- 
ing femtosecond pulses. 

There are still many open questions to answer. In 
this work, only the pulse shaping in the time domain is 
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considered. The spatial resolution of CARS microscopy 
is diffraction limited, without employing quantum 
effects. Potma et al. apply the concept of focus 
engineering^^^^ to enhance the sensitivity of CARS 
microscopy to "chemical interfaces". However, coherent 
control strategies by pulse-shaping in the spatial domain 
to improve the spatial resolution of CARS microscopy 
have not been given much attention. Another point is 
chemical selectivity in CARS microscopy. Experiments 
have demonstrated that not only can the CARS signals 
from different vibrational modes in one molecule inter- 
fere but also from different molecules. By adjusting the 
pulse phases, the signal from one special molecule can 
be effectively enhanced while that of another molecule is 
suppressed^^. Thus, the simultaneous enhanced chemical 
selectivity and spatial resolution of CARS microscopy 
needs to be explored systematically in the future by 
phase shaping in the time and spatial domains. 



The variational method is used to analyze all the sta- 
tionary points for |P^p with respect to the phase function 
^pr{ujpr)' The variational condition S \Prf = requires 



From this criterion, we have 



(A5) 



Re[ 



exp \i (<l>p^(cjp^) + a(cjp^))]] = 



(A6) 



which is equivalent to 
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Appendix A: Variation of resonant Signal 

In the appendices, the following notations are used for 
clarity 



sin {^pr{uj) + = a2{cjo) 

cos {^pr{uj) + a{uj)) = 62 (cj) 
A2 



1 



_ 3x^ 

e 2A^ 



_ 3x'^ 

e 2A^ 



sin {^pr{x) -\- a {x)) dx 



Vx^ + r2 
1 

Vx^ + r2 



cos {^pr{x) + a{x)) dx>=B2 



and 



{/: 



sin($pr(w)) 


= ai(w) 


(Al) 


cos ($pr(w)) 


= &iM 


(A2) 




= Ai 


(A3) 


2a7 cos {^pr{x)) dx > 


= Bi 


(A4) 



SO the general extremal phase functions ^pr{ojpr) for all 
the critical points of |PrP satisfy 



^pr{ujpr) = L[ujpr]7r + arctan(c(jp^/r) + constant^ 

where L[ujpr] G {0, 1} is the integer function of ujpr. It is 
easy to find that when L[ujpr] = or 1, i.e. ^pr{oJpr) = 
arctan(a;pr/r) + constant, then |P^p reaches its global 
maximal value. When L[ijjpr] = IIeaviside(cc;pr), i.e. 
^pr{oJpr) = arctan(u;p^/r)+Heaviside(cc;p^)-7r-hcon5tant, 
|P^P has its minimal value 0. Hence, the arctan(co'p^/r) 
phase function is the optimal solution for maximizing 
the peak resonant signal intensity. For the other cases 
of K\ (jjpr]^ the phase functions are either local extremal 
points or saddle points. However, numerical optimization 
did not reveal any local extremal points, which may indi- 
cate that the control landscape^^ for the resonant signal 
intensity is trap-free in this pure phase shaping strategy. 
A detailed proof of the trap-free landscape for CARS sig- 
nal needs to be investigated. 



Appendix B: Variation of the difference between the 
Signal and Background 

For simplicity in the proof, we set P^o = Pr/C 
and Pnro = Pnr/Xnr^ and introduce a parameter A = 
^ (Xnr/C')^, then the optimization of the objective func- 
tional 



J — |Pro|^ — A |Pnro|^ 



(Bl) 



is equivalent to optimizing the original functional Eq. 
(14), because \Prf - k |P„,|' = C2(|P,o|' - A |P„.o|'). 
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Following a similar procedure to optimizing |PrP in 
Appendix A, the optimal phase function for J satisfies 



\{aiBi - Aihi) = a2B2 - 



(B2) 



hence 



B2= fZa^ ^ ^^^^{x-COs{^pr^{x))-r-sin{^pr^{x)))dx^((B7) 



So we have 



tan {^pr{cjo)) = 



bi{uj) A(a;2 + r2)5i - B2UJ - A2T 

(B3) 

According to the expressions for \Pr^\ojas)\'^ and 
Pnr\oJas)\'^ ^ the valuc of J will not change if a trivial 
phase constant is added to ^pr{ojpr). So the general so- 
lution for maximizing J can be expressed as the sum of 
one special solution and a trivial phase constant. The 
numerical simulation in Fig. (6) shows that the optimal 
phase is just an odd function of ujpr. So ^pr{ujpr) to opti- 
mize J just makes sin (<l>pr(^pr)) an odd function of uOpr 
and cos (^pr(^pr)) an even function of ujpr , thus 



Ai= e 2A^ sin{^pr{x))dx^O 



(B4) 



(B5) 



^2= /r^e i&"^^^(cc-sin(<l>p,(cc))+rcos(<^>p,(x)))da;/0(B6) 



The special solution ^p^iuj) = arctan( ^ — 

satisfies the assumption, so the general solution for max- 
imal J is 



^pr(^) = arctan( 



r-Af^(c.2 + P) 



) + constant 



Here we define 7 = ^1/742, which could not be de- 
termined analytically and was determined by iteratively 
solving the equation 



roo 2A2 cos( arctan T ^ ^^^dcci 



— , ^ sin( arctan ( ^ ]+cx(x) ]dx > 



(B9) 

Given parameters A ,A and F, we can always get the 
numerical value of 7. 
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